home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / public / pbmplus / pgm / pgmhist.c < prev    next >
C/C++ Source or Header  |  1994-08-01  |  2KB  |  86 lines
  1. /* pgmhist.c - print a histogram of the values in a portable graymap
  2. **
  3. ** Copyright (C) 1989 by Jef Poskanzer.
  4. **
  5. ** Permission to use, copy, modify, and distribute this software and its
  6. ** documentation for any purpose and without fee is hereby granted, provided
  7. ** that the above copyright notice appear in all copies and that both that
  8. ** copyright notice and this permission notice appear in supporting
  9. ** documentation.  This software is provided "as is" without express or
  10. ** implied warranty.
  11. */
  12.  
  13. #include "pgm.h"
  14.  
  15. void
  16. main( argc, argv )
  17. int argc;
  18. char *argv[];
  19.     {
  20.     FILE *ifp;
  21.     gray maxval, *grayrow;
  22.     register gray *gP;
  23.     int argn, rows, cols, format, row;
  24.     int i, *hist, *rcount, count, size;
  25.     register int col;
  26.     char *usage = "[pgmfile]";
  27.  
  28.     pgm_init( &argc, argv );
  29.  
  30.     argn = 1;
  31.  
  32.     if ( argn < argc )
  33.     {
  34.     ifp = pm_openr( argv[argn] );
  35.     argn++;
  36.     }
  37.     else
  38.     ifp = stdin;
  39.  
  40.     if ( argn != argc )
  41.     pm_usage( usage );
  42.  
  43.     pgm_readpgminit( ifp, &cols, &rows, &maxval, &format );
  44.     grayrow = pgm_allocrow( cols );
  45.  
  46.     /* Build histogram. */
  47.     hist = (int *) malloc( ( maxval + 1 ) * sizeof(int) );
  48.     rcount = (int *) malloc( ( maxval + 1 ) * sizeof(int) );
  49.     if ( hist == (int *) 0 || rcount == (int *) 0 )
  50.     pm_error( "out of memory" );
  51.     for ( i = 0; i <= maxval; i++ )
  52.     hist[i] = 0;
  53.     for ( row = 0; row < rows; row++ )
  54.     {
  55.     pgm_readpgmrow( ifp, grayrow, cols, maxval, format );
  56.         for ( col = 0, gP = grayrow; col < cols; col++, gP++ )
  57.         hist[(int) *gP]++;
  58.     }
  59.  
  60.     pm_close( ifp );
  61.  
  62.     /* Compute count-down */
  63.     count = 0;
  64.     for ( i = maxval; i >= 0; i-- )
  65.     {
  66.     count += hist[i];
  67.     rcount[i] = count;
  68.     }
  69.  
  70.     /* And print it. */
  71.     printf( "value\tcount\tb%%\tw%%\n" );
  72.     printf( "-----\t-----\t--\t--\n" );
  73.     count = 0;
  74.     size = rows * cols;
  75.     for ( i = 0; i <= maxval; i++ )
  76.     if ( hist[i] > 0 )
  77.         {
  78.         count += hist[i];
  79.         printf(
  80.         "%d\t%d\t%5.3g%%\t%5.3g%%\n", i, hist[i],
  81.         (float) count * 100.0 / size, (float) rcount[i] * 100.0 / size );
  82.         }
  83.  
  84.     exit( 0 );
  85.     }
  86.